Stripes in Doped Antiferromagnets: Single-Particle Spectral Weight 
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Recent photoemission (ARPES) experiments on cuprate superconductors provide important guide- 
lines for a theory of electronic excitations in the stripe phase. Using a cluster perturbation theory, 
where short-distance effects are accounted for by exact cluster diagonalization and long-distance ef- 
fects by perturbation (in the hopping), we calculate the single-particle Green's function for a striped 
t-J model. The data obtained quantitatively reproduce salient (ARPES-) features and may serve 
to rule out "bond-centered" in favor of "site-centered" stripes. 

PACS numbers: 74.72.-h, 79.60.-i, 71.27.+a 



There is by now substantial experimental evidence 
for a tendency of doped holes in high-temperature su- 
perconductors (HTSC) to form stripes, leaving behind 
locally antiferromagnetic (AF) domains. The stripes can 
be static like in Lai.48Ndo.4Sro.i2Cu0 4 (Nd-LSCO) or 
dynamic like in YBa2Cu307_ x [[|. At present, there is 
an intensive discussion and controversy, whether stripes 
are directly connected with and beneficial for the mi- 
croscopic mechanism in HTSC [^|— pTT| . What is clear, 
however, is that the apparent presence of static or dy- 
namic stripes crucially influences low-energy excitations 
and thus the foundations for such a microscopic theory. 
Evidence for this has recently been accumulated by angle- 
resolved photoemission spectroscopy (ARPES) both on 
the static stripes in the Nd-LSCO system || and on dy- 
namic domain walls in the LSCO compound |l2|]l3| ]. The 
electronic structure revealed by ARPES contains char- 
acteristic features consistent with other cuprates, such 
as the flat band at low energy near the Brillouin zone 
face (k = (tt, 0)). In Nd-LSCO, the frequency-integrated 
spectral weight is confined inside ID-segments in k-space, 
deviating strongly from the more rounded Fermi surface 
expected from band calculations. 

In this Letter, we present a numerical study of the single- 
electron excitations in a striped phase via cluster per- 
turbation theory (CPT) and a detailed comparison 
with recent ARPES results. The basic idea of our appli- 
cation of the CPT is indicated in Fig.(|l|): it is based on 
dividing the 2D plane into alternating clusters of metal- 
lic stripes and AF domains. The individial clusters are 
modeled by the microscopic t-J hamiltonian and solved 
exactly via exact diagonalization (ED). Then the inter- 
cluster hopping linking the alternating metallic and AF 
domains is incorporated perturbatively via CPT on the 
basis of the exact cluster Green's functions thus yield- 
ing the spectral function of the infinite 2D plane in a 
striped phase. Using this CPT approach, the important 
short-distance interaction effects within the stripes are 
taken into account exactly while longer-ranged hopping 
effects are treated perturbatively. A study of the prop- 
erties of experimentally observed stripe phases solely by 
ED Q is precluded by the prohibitively large unit cells. 



The manageable clusters for ED are simply too small to 
accommodate even a single such unit cell. 

Our main results are: (i) close to k = (7r, 0) we see, like 
in experiments, a two-component electronic feature (see 
Fig.(^)): a sharp low-energy feature close to Ep and a 
more broad feature at higher binding energies. Both fea- 
tures can be explained by the mixing of metallic and an- 
tiferromagnetic bands at this k-point. (ii) the excitation 
near (7r/2, tt/2) is at higher binding energies than the low- 
energy excitation at (7r, 0) and of reduced weight, (hi) 
the integrated spectral weight of the cluster-stripe cal- 
culation resembles the quasi-one-dimensional segments 
in momentum space (see Fig. (Ha)) as seen in the Nd- 
LSCO experiment. Also in agreement with the Nd-LSCO 
experiment our calculation finds the low-energy excita- 
tions near (±7r, 0) and (0, ±7r) (Fig.([|b)). Interestingly, 
this agreement with experiment occurs only for so called 
"site-centered" metallic stripes (as shown in Fig.(Q)) and 
not for "bond-centered" metallic stripes. This seems im- 
portant since it has been argued that, for bond- 
centered stripes, superconductivity is expected to survive 
stripe ordering. In the DMRG calculations by White 
and Scalapino |l7j as well as in dynamical mean-field 
(DMFT) studies by Fleck et al @, bond-centered and 
site-centered domains are very close in energy (ground- 
state). However, our technique allows to distinguish 
them dynamically. 

The computational technique for our calculation of the 
single-particle spectral weight A(k, uS) and the Green's 
function is illustrated in Fig. ([]]). We solve the AF clus- 
ter (here, N x 3) and the metallic cluster (here N x 1 
with a hole filling rih = 0.5) by ED and combine the in- 
dividual clusters to an infinite lattice via CPT as in the 
homogeneous case. Where it was technically feasible, we 
extended the unit cell by a factor of two and diagonal- 
ized two AF 3-leg clusters with a staggered magnetic field 
pointing in opposite directions, resulting in a 7r-phase 
shifted Neel order in the final configuration. This site- 
centered "3 + 1" configuration with 7r-phase shifted Neel 
order of stripes was first suggested by Tranquada et al. 
[Q. Bond-centered stripes, on the other hand, are mod- 
eled by 2-leg ladders with alternating filling (half-filled, 
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rih = and doped, = 1/4). In the following we will 
refer to this bond-centered configuration as " 2 + 2" . 
Holes can propagate out of the metallic stripes into the 
AF insulating domains via the inter-cluster hoppings. In 
the Hubbard model, these hoppings correspond to one- 
body operators and can be treated within a systematic 
strong-coupling perturbation expansion. For homoge- 
neous systems such an expansion was constructed in ref. 
[]l4]l. We take the lowest order contribution: 
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which is of RPA form. Here, P is a super lattice wave 
vector and is the Green's function of the "oo-size" 
2D system, however, still in a hybrid representation: 
real space within a cluster and Fourier-space between 
the clusters. This is related to the fact that Goo(P,z) 
is now an M x M matrix in the space of site indices 
(in the inhomogeneous stripe configuration of Fig.(|l|) 
M = 7Vx3 + 7Vxl = AN). Likewise, e(P) and G duster 
are M x M matrices in real space with e(P) standing 
for the perturbation, which includes hoppings out of the 
clusters. A true Fourier representation of Goo in terms 
of the original reciprocal lattice then yields the lowest 
order CPT approximation. As discussed in ref. 
the homogeneous case, the approximation in eq.(| 
act for vanishing interaction. When the interactions are 
turned on, eq.(|l|) is no longer exact (apart from the lo- 
cal limit, i.e. t = t' = 0), but strong interactions are 
known to be important mainly for short-range correla- 
tions. These correlations are incorporated with good ac- 
curacy in modest-size clusters and are treated here by ED 
within the cluster. It has been shown by Senechal et al. 
Q that the CPT reproduces the spectral weight of the 
ID and 2D Hubbard models in quantitative agreement 
with exact results. 

To allow for larger cluster sizes A/", we diagonalize the t-J 
model and take its spectral (Green's) function as our lo- 
cal G c i us ter in ec i-(|l|) as an approximation to the Hubbard 
model's (cluster-) Green's function. A comparison of the 
Hubbard and t-J model's spectral function on small clus- 
ters [[l9j shows that the strong low energy peaks have 
similar dispersion and weight in both models, the main 
difference being a transfer of incoherent high energy spec- 
tral weight from momenta near (7r,7r) to (0,0). The t-J 
hamiltonian is defined as 
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The hopping matrix element tij is nonzero only for near- 
est (t) and next-nearest neighbors (£'). The second sum 
counting the Heisenberg interaction J runs over all near- 
est neighbor pairs. No double occupancy is allowed. 
We have chosen commonly accepted values for the ra- 
tio t'/t = -0.2, J jt = 0.4, where t « 0.5eV. In this 



Letter we present calculations for systems with N = 8 
based on diagonalizations of a N x 3 = 24-site half-filled 
3-leg ladder and a quarter-filled 8-site chain. Results for 
smaller N = 6 do not differ much from N = 8 results. 
N = 6, however, is somewhat pathological, since it has 
an odd number of electrons in the quarter-filled chain. 
We proceed to the discussion of the spectra: Fig.^a) 
shows the experimental ARPES results for LSCO at the 
superconductor-insulator transition (doping x = 0.05) 
]l2[ . To enhance the structure in the obtained spectra, 
the authors of ref. jl2| plotted the second derivative of 
the ARPES spectrum, so areas with high second deriva- 
tive are marked white and areas with low curvature (i.e. 
flat intensity) are black. This result is compared with 
the theoretical CPT calculation for different stripe con- 
figurations with overall doping of x = 1/8. Figs.(H|b,c) 
are for the "3 + 1" site-centered configuration. Fig.(||b) 
shows the result for a "3 + 1" configuration with alter- 
nating (i.e. 7r-phase shifted) Neel order between the 3-leg 
ladders (induced by a staggered magnetic field B = O.lt) 
without next-nearest neighbor hopping, Fig.(||c) shows 
the result for the "3 + 1" stripe configuration with next 
nearest-neighbor hopping t' = — 0.2t, however without 
Neel order (due to the reduced symmetry, a t' diagonal- 
ization with staggered field is not technically feasible). 
Fig.(||d) shows the result for a bond-centered "2 + 2" 
stripe configuration. We observe that the spectra for the 
site-centered "3 + 1" configuration (Fig.(||b,c)) are in sur- 
prisingly good agreement with experiment, and, that the 
t' — calculation (Fig.(|]b)) results in much more co- 
herent bands due to the enforced Neel order in the 3-leg 
ladders. Similar to recent DMFT calculations |18|, the 
sharp excitation near the Fermi surface around (n, 0), 
that has been interpreted by Ino et al. as the quasiparti- 
cle peak in the SC state, is visible as well as a dispersive 
band at higher binding energies which (at least away from 
(7r, 0); see discussion below) can be interpreted as rem- 
nants of the insulating valence band resulting from the 
AF domains. Especially in the Neel ordered configura- 
tion, we observe a very coherent and pronounced band. 
This clear dispersion is also visible in the (7r, tt) direction 
near (7r/2,7r/2). We note that, in agreement with the 
experimental result (Fig.(^a)), the excitation at (jr, 0) is 
at significantly lower binding energy than the excitation 
at (7r/2,7r/2). Neither the 2D t-J model nor a 2D t-t'- 
t"-J model, with its parameters fitted to the insulating 
state, can reproduce this result [2C| ]. This is a crucial ef- 
fect of the stripe assumption: With the stripes oriented 
along the ^-direction, the metallic band is dispersionless 
in x-direction. Therefore, at k = (71" , 0), the minimum of 
the metallic spinon band (located at (k x ,0) for any k x ) 
hybridizes with the top of the insulating valence band 
resulting in a two-peak structure with one peak pushed 
to higher and the other pushed to lower binding energies 
(see below). At k = (7r/2,7r/2), on the other hand, the 
metallic band has crossed the Fermi surface (its hp being 
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7r/4) and no mixing takes place. Finally, the "2+2" bond- 
centered stripe configuration (Fig.(||d)) does not show 
much resemblance to the experimental result. Its main 
band is much more two-dimensional, normal metal-like, 
comparable to the dispersion of a 2D tight-binding band. 
Fig.(||a) plots the integrated spectral weight n(k) for the 
"3 + 1" site-centered stripe configuration with t' = 0. 
Although not as clear as in the Nd-LSCO ARPES ex- 
periment (from ref. ||), the "Fermi surface" is rather 
one-dimensional in structure. Like in Nd-LSCO the low 
energy excitations (shown in Fig.(|3|b), calculated by in- 
tegrating over a Auj = 0.2t window below the Fermi- 
energy for each k-point in the Brillouin zone) are located 
near the (±7r,0), (0, ±tt) points in momentum space. In 
Fig.(||b) we notice the 8x8 square lattice of bright points. 
This is the repeated Brillouin zone of the supercell con- 
sisting of(3 + l + 3 + l)x7V = 8x8 lattice sites (due to the 
Neel order in x-direction). Clearly, the low energy exci- 
tations in the momentum space of the supercell are near 
the (±7r, 0) and (0, ±tt) points as well. Fig.(||c) shows the 
integrated spectral weight n(k) of the "2 + 2" t-J stripe 
calculation. Here, the CPT "Fermi surface" is much more 
rounded, similar to the quasi-2D Fermi surface known 
from band calculations. The low energy excitations are 
located isotropically around the "Fermi surface" as well 
(Fig.(|d)). The loss of one-dimensionality observed for 
this bond-centered " 2 + 2" stripe configuration is accom- 
panied by a substantial enhancement of spectral weight 
near (7r/2,7r/2). From these A(k, uS) results we conclude 
that, at least in the Nd-LSCO system, the stripes are 
site-centered and of "3 + 1" type. 

With our technique we are able to resolve for each exci- 
tation, whether its main origin is from the insulating or 
the metallic part of the stripe configuration: In Fig.(Q) 
we compare the spectra for the unperturbed stripe con- 
figuration with inter-cluster hopping set equal to zero (in 
Fig.(|]a): solid curve for the AF domains, shaded curve 
for the ID metal) with the result of our CPT calculation 
(solid line in Fig.([I]b)) with inter-cluster hopping t (for 
the Neel ordered "3 + 1" configuration). All of the spec- 
tral weight in inverse photoemission (u > 0) naturally 
stems from the chain since the half- filled 3-leg t-J ladder 
does not have target states for an inverse photoemission 
process. The stripes are oriented along the ^-direction. 
Therefore, we can conclude that peaks (in the calcula- 
tion with t = 0), that show a dispersion along (0,0) to 
(7r, 0) direction stem from the AF domains (solid line in 
whereas the metallic excitations prior to the mixing 
(shaded curve in (f|a)) are dispersionless. By comparing 
the three curves, we therefore conclude that the sharp 
quasiparticle peak near (7r, 0) results from the mixing 
of the (dispersionless) (k x ,0) minimum of the metallic 
spinon band and the top of the insulating valence band 
situated at (7r,0). Going from (-zr, 0) to (7r,7r), the metal- 
lic band becomes dispersive and crosses, in agreement 
with experiment [pLH the Fermi surface at k y — ir/A 



(since it is quarter-filled). The dispersion of the insu- 
lating band, however, is in the opposite direction. For 
this reason, in the final spectrum, we observe that the 
sharp quasiparticle peak at (7r, 0) becomes dispersive go- 
ing into (0, 7r) direction and eventually crosses the Fermi 
surface, however, with diminishing weight due to the ab- 
sence of mixing with the insulating band. This effect is 
best visible in the grayscale plot of Fig.([|b). This finding 
may serve to clarify questions raised in the experimental 
ref. m concerning the origin of the quasiparticle peak at 

(M). 

To summarize, the single-particle spectral weight A(k, uo) 
was calculated for different stripe configurations employ- 
ing an application of the cluster perturbation technique 
for inhomogeneous systems. This technique allows to ob- 
tain all k-points and, therefore, allows for a detailed com- 
parison with ARPES data. The A(k, uo) results for the 
"3 + 1" site-centered configuration display salient features 
observed in experiments such as a two-peak structure 
around (7r, 0) with a sharp excitation close to the Fermi- 
energy and a broader feature at higher binding energies, 
a quasi- ID distribution of spectral weight, and low en- 
ergy excitations located around the (7r, 0)-points in the 
Brillouin zone. The theoretical results suggest to rule out 
the alternative bond-centered " 2 + 2" configuration. 
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(05SB8WWA1) and DFG (HA1537/17-1). The calcula- 
tions were carried out at the high-performance comput- 
ing centers HLRS (Stuttgart) and LRZ (Mimchen). 



[1] J.M. Tranquada et a/., Nature 375, 561 (1995). 
[2] H.A. Mook et a/., Nature 404, 729 (2000); Nature 395, 
580 (1998). 

[3] X.J. Zhou et al, Science 286, 268 (1999). 

[4] R.P. Sharma et a/., Nature 404, 736 (2000). 

[5] N. Ichikawa et a/., cond- mat /99 10037 (1999). 

[6] V.J. Emery, S.A. Kivelson, O. Zachar, Phys. Rev. B 56, 
6120 (1997); S.A. Kivelson, E. Fradkin, V.J. Emery, Na- 
ture 393, 550 (1998), and ref. therein. 

[7] J. Zaanen, J. Phys. Chem. Sol. 59, 1769 (1995). 

[8] C. Castellani, C. Di Castro, M. Grilli, Z. Phys. B 103, 
137 (1997). 

[9] M. Vojta, S. Sachdev, Phys. Rev. Lett. 83, 3916 (1999). 
[10] S.R. White, D.J. Scalapino, Phys. Rev. B 60, R753 
(1999). 

[11] P. Wrobel and R. Eder, Phys. Rev. B 62, 4048 (2000). 
[12] A. Ino et a/., Phys. Rev. B 62, 4137 (2000). 
[13] A. Bianconi et a/., Phys. Rev. Lett. 76, 3412 (1996). 
[14] D. Senechal et al, Phys. Rev. Lett. 84, 522 (2000). 
[15] T. Tohyama et a/., Phys. Rev. Lett. 82, 4910 (1999). 
[16] Y.A. Krotov, D.-H. Lee, A.V. Balatsky, Phys. Rev. B 56, 
8367 (1997). 

[17] S.R. White, D.J. Scalapino, Phys. Rev. Lett. 80, 1272 



3 



(1998); Phys. Rev. Lett. 81, 3227 (1998); Phys. Rev. B 

55, 6504 (1997). 
[18] M. Fleck et al, Phys. Rev. Lett. 84, 4962 (2000). 
[19] H. Eskes, R. Eder, Phys. Rev. B 54, R14226 (1996). 
[20] A. Nazarenko et al, Phys. Rev. B 51, 8676 (1995); V.I. 

Belinicher et a/., Phys. Rev. B 54, 14914 (1996); R. Eder 

et a/., Phys. Rev. B 55, R3414 (1997); F. Lema and A. 

Aligia, Phys. Rev. B 55, 14092 (1997). 



i i i-.-i.-i i i 



N 



T 



T 



FIG. 1. Visualization of the cluster perturbation approach 
for stripes: the groundstates for the half-filled three-leg lad- 
der (3 x N) and the quarter- filled 1-leg chain (1 x N) are 
calculated exactly via exact diagonalization. The alternating 
clusters are then coupled via the inter-cluster hopping which 
is treated perturbatively. 




FIG. 3. Integrated spectral weight of site-centered (a,b) 
and bond-centered (c,d) stripe configurations; (a,c) total in- 
tegrated weight in photoemission (n(k)), (b,d) low energy ex- 
citations (integrated weight in Ef — 0.2t < uj < Ef)- The 
data are plotted for the whole Brillouin zone with the T-point 
in the center. The result of the stripe calculations have been 
symmetrized to account for the differently oriented stripe do- 
mains in real materials. Regions of high spectral weight cor- 
respond to white areas. 
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FIG. 2. ARPES results for La 2 -xSr x Cu04, theoretical sin- 
gle particle spectral function A(k, cj): (a) displays experimen- 
tal ARPES results by Ino et al. |L2]. The grayscale cor- 
responds to the second derivative of the original measured 
data. Flat regions are black, regions with high curvature (i.e. 
peaks) are white. (b),(c),(d) show the results of the stripe 
CPT calculation for A(k,uj) (b: "3 + 1" site-centered \t' — 0, 
c: "3 + l";t' = -0.2t, d: "2 + 2" bond-centered;*' = 0). 
Here, A(k,uj) is plotted directly with maximum intensity cor- 
responding to black. 
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FIG. 4. Single particle spectral weight prior and after ap- 
plication of the cluster perturbation theory: In (a) the shaded 
curve gives the spectrum of the ID metallic chains and the 
solid line corresponds to the spectrum of the uncoupled AF 
domains (see Fig.Q). In (b) the result of the cluster perturba- 
tion is plotted. 



